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Abstract. - We study the upper critical fleld of the A and B phases in the triplet supercon- 
ductor PrOs4Sbi2 within the p+h-wave superconductivity proposed recently for this material. 
The present result is compared with Hc2(t) and H*(t), the boundary between the A and B 
phase in PrOs4Sbi2, reported earlier and with more recent data of Hc2(i) for the single phase 
crystal. We find Hc2(t)'s for both the two phase crystal and the single phase crystal are de- 
scribed by the model for the A phase. From this fltting one can deduce the Fermi velocity as 
V = 2.5 X 10® cm/s. On the other hand Hc2(t) for the B phase is found to be somewhat smaller 
than H*(t), which is rather puzzling. 



Introduction. - Superconductivity in the filled skutterudite PrOs4Sbi2 discovered in 2002 
by Bauer et al [1-3] has generated a big sensation. First the presence of at least two distinct 
phases (A and B phases) in a magnetic field, with both gap functions having point nodes 
is surprising. Further these superconductors belong to the triplet pairing with broken chiral 
symmetry [4-6]. However, the exact location of the A-B phase boundary is still controversial 
[7]. More recently Measson et al [8,9] have discovered many PrOs4Sbi2 crystals with only a 
single phase. Therefore the nature of this single phase has to be addressed. In this paper we 
assume that the gap functions of the A and B phases of PrOs4Sbi2 are given as [6] 



A^(k) ^ dAe±'^--(l ~kt~ kl - kt)). (1) 
AB(k) = dBe±*^Hl-fc') (2) 



Here e^*"^' is one of e"^^ = {ky + ik,)/^kl + kl,e">'^ = {k, + ik^)/^^k^^ + kl,e''^^ = (k^ + 

iky)/ \Jk'i + ky. The factor of 3/2 ensures proper normalization of the angular dependence 
of the order parameter. In Eq.(2) the nodal direction is chosen to be parallel to (GDI). We 
note that ]Ayi(k)] retains the cubic symmetry while ]As(k)] has only the axial symmetry. 
In the absence of an external perturbation, we assume that the nodal direction of ]AB(k)| is 
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parallel to H, the magnetic field, when H is parallel to one of the crystal axes. Although many 
other gap functions for PrOs4Sbi2 have been proposed so far, none can describe the thermal 
conductivity data of Izawa et al as discussed in [4,6]. For a discussion of our methodology for 
determining the gap functions of nodal superconductors please see Ref. [10]. 

In the following we study the upper critical field Hc2(T) of the A and B phases of PrOs4Sbi2 
in terms of the linearized gap equation [11]. For imconventional superconductors Gor'kov's 
original formulation must be modified as shown by Luk'yanchuck and Minecv [12]. The upper 
critical field thus obtained for the A phase describes consistently Hc2 of the A phase as obtained 
by Measson et al [7]. On the other hand H*(T) is somewhat larger than Hc2(T) in the B phase 
at low temperatures. 

Of course H*(T) is not exactly equal to Hc2(T) of the B phase. But it is rather puzzling 

that H*(T) is larger than the predicted Hc2(T) of the B phase. Meanwhile, some of the 
single crystals exhibit only a single phase transition. Measson has reported Hc2(T) of the 
single phase system at SCES05 in Vienna [8]. The observed Hc2(T) is practically the same as 
H*(T), and consistent with the one for the A phase. This, perhaps, implies that the A phase 
is more robust than the B phase. Also, perhaps H*(T) is associated with the A phase. 

Prom the spatial configuration of A(r,k) at H=Hc2(T), we can deduce the stable vortex 
lattice structure as in [13, 14]. Except in the immediate vicinity of the superconducting 
transitions (T < 0.95Tc) we find the square vortex lattice is more stable than the hexagonal 
one for H || [001] or any of the cubic axes, while the hexagonal lattice is more stable for 
H )| [111] for example. Recently the vortex lattice in PrOs4Sbi2 was studied via small angle 
neutron scattering (SANS) by Huxley et al [15]. They find only a distorted hexagonal lattice 
at T=100 mK and H =0.1 T, most likely due to poor sample quality. 

Upper critical field. - As in the chiral p-wave superconductor the upper critical field is 
determined from the linearized gap equation for A(r, k): 

A(r,k)~ (e"^ + Ce-'*(at)2)lO> (3) 
where [0 > is the Abrikosov state [16]: 

|0> = ^c„e-^-^^'-"'=(^+'2')-^ (4) 

n 

and = --j=^{—idx — dy + 2ieBx) is the raising operator and we have assumed H || c. 
Then following Refs. [13, 14] the upper critical field is determined by 

- Ini = ^(1- < f (/^e-^»^(i-^) X (1 + 2Cpu^l - z^)))) (5) 

-C\nt = r ^(1- < f {fe-P-'i^--') X {pu\l - z^) + 
Jo smh u 

C(l - Apu\l - z") + 2p^u\\ - z'f)))) (6) 

where t= ^,p= ^^^yp^, where v is the Fermi velocity and 

/ = l-cos''6'-sin''6'(sin'*0 + cos''(?!)), (A - phase) (7) 
= l-cos*6», (B- phase) (8) 



D. Parker, K. Maki and H. Won: Upper critical field Hc2 in PROS4SB12 



3 



and (. . .) means / df2/(47r) — Also < >= ^ and || for the A and B phases respectively. 



45 

Then for the A phase we obtain the following asymptotics for f ~ 1: 

-Ini = ^^p(l-2C) and (9) 
-Clnt = ^^p{-l + 5C) (10) 

where C(3) = 1.202.... From this wc find C = - 1 and p = 7^(1 + ^Jl){-lnt) = 
0.5942(- \nt). On the other hand, for t=0 we find 

= ln(47po)+ < f if ln(l - z^)) - 2C and (11) 
= C(ln(47po)+ < f if ln(l -z^)+2)-l (12) 

where 7 = 1.780... is the Euler constant and < f >-^< - z'^) >= -0.40104. These 

yield C= i(V3- 1) and 

From this we obtain 

MO) = = 0-80696 (14) 
dt l*=i 

Hc2(t) and C{t) arc evaluated numerically and shown in Fig. 1 and Fig. 2 respectively. Also 
in Fig. 1 Hc2(t) from the A phase is compared with the data from Ref. [7]. As is readily seen, 
wc have an excellent fit, where we used v = 2 Ah x 10^ cm/sec which is very reasonable. This 
V is somewhat smaller than the one deduced from the thermal conductivity data [6]. However, 
it is known that the thermal conductivity data in the clean limit (and not in the superclean 
limit) is less sensitive to the actual value of v. In the same way v is smaller than that deduced 
from de Haas- van Alphcn data [17]. But in heavy- fermion systems such as UPta, it is well 
known that dHvA gives a lighter effective mass than that deduced via the specific heat. 
For the B phase we find for t ~ 1 

-\nt = ^i^p(l - 2C) and (15) 
-Clnt = 5lMp(_i + 5C) (16) 

which give C = - 1 and p = 31^(1 + ^Jl){- Ini) 0.5362(- Ini). For t ~ 0, on the 
other hand, we have the same set of equations as in Eqs.(lO) and (11), where 

< f >"^< f ln(l - z"^) >= -0.25973 (17) 

This gives C =\{V^-l) and 

po = exp[1.2295025] = 0.3786 (18) 
47 
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We then find 



Both Hc2(i) and C(t) for the B phase are obtained numerically and shown in Fig. 1 and Fig. 
2 respectively as before. Also H*(t) from Ref. [7] is shown in Fig. 1 as well, though there is no 
reason that H*(t) should correspond to Hc2(t) of the B phase. Indeed we see that Hc2(t) for 
the B phase is substantially smaller than the observed H*(t), especially at low temperatures. 

We also compare in Fig. 1 the observed Hc2(t) of the single phase crystal with previous 
data, as well as our predictions for the A and B phases. The single phase data appears 
essentially identical to H*(t) of the two phase data, suggesting that the single-phase crystals 
have only the A phase. H*(t) and Hc2(t) of the single phase crystal are fully consistent with 
Hc2(t) of the A phase, if Tc is taken as 1.72 K, as indicated by the middle fit line in Fig. 1. 
Such a reduction in Tc is readily reproduced in the presence of impurity scattering [12]. Also 
we recall the scaling of Hc2(T/Tc, r)/Hc2(0, F) in the presence of impurity scattering. These 
strongly suggest that both H*(t) and Hc2(t) of the single phase system should be associated 
with the A phase. As we have mentioned earlier the actual A-B phase boundary is still 
controversial; below in Fig. 3 we show a measurement [4,6] indicating a phase boundary lying 
much lower in the H-T plane, and hence more consistent with our preduction for Hc2(i) of the 
B phase. Clearly experiments with cleaner crystals are desirable. 

Also from Fig. 2 we notice that C(t) for the B phase is practically the same as the one 
obtained for the A phase. Further the present C(t) is practically the same as C(t) obtained 
for the chiral p-wave superconductor [13]. We believe there should be an analytical reason 
that C(t) is universal for the class of chiral superconductors. This means that the Abrikosov 
state consists of the n=0 Landau wave function with an admixture of the n=2 Landau wave 
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Fig. 2 - C(T) for the A and B phases. 
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Fig. 3 - Alternate phase diagram of PrOs4Sbi2 



functions with opposite chirality. Further the coefRcient of the n=2 Landau wave function is 
a universal function of < = T/Tc- 

Vortex lattice. - A relatively large C (i.e. C > — 1 ~ 0.22474) appears to be 
characteristic of triplet superconductivity with chirality ±1, as observed in Sr2Ru04. This 
suggests strongly that the square vortex lattice is more stable than the usual hexagonal lattice 
except in the immediate vicinity of the superconducting transitions at T^a and T^b- 
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In order to examine the vortex lattice structure we calculate the Abrikosov parameter [13, 14] 

_ (|A(r)n 

For both the A and B phases in a magnetic field H || c, we obtain 

I3a = J^ \^(_l)"me-^^("'+™') X (1 + 2C)-2 X / dxe'^^'il + C^{f{x + x^) + f{x - xi)) x 

{fix + X2) + fix - X2)) + C^fix + Xi)fix - Xi)fix + X2)fix - X2m) 

where fix) = x"^ — 1,R = tan(^), xi = \J ■^in — m), X2 = \J^in + m) and is the 

apex angle of the vortex lattice. Here R = ^/2> corresponds to the regular hexagonal lattice, 
while i? = 1 is the square lattice. The above expression is the same as in the chiral p-wave 
superconductor discussed in [13]. As is well known, for the most stable lattice (3a takes the 
minimum value. From the analysis of the chiral p-wavc superconductor [13] we conclude that 
the square vortex lattice is the most stable for T < 0.95Tca and T < 0.95Tcb in the A and B 
phase respectively. 

Of course the present analysis applies only in the vicinity oi H = Hc2{T). In order 
to find the phase boundary between the square vortex and the hexagonal vortex lattice a 
parallel analysis, as was done for the vortex lattice for d-wavc superconductivity by Shiraishi 
et al [18], is required. But we expect the transition line should be around k^^Hc2{T) where 
K is the Ginzburg-Landau parameter of PrOs4Sbi2. Making use of the magnetic penetration 
depth determined by Chia et al [19] we can extract k ~ 30.0. This means that the square 
vortex lattice should be visible for H= 0.1 T in an ideal crystal. As already mentioned Huxley 
et al [15] studied the vortex lattice via small angle neutron scattering (SANS), but saw only a 
distorted hexagonal lattice. In a magnetic field H ]] [001], most of the nodal excitations in both 
the A and B phases of PrOs4Sbi2 are moving parallel to the [100] and [010] directions. Hence 
a square vortex lattice aligned parallel to the crystal axis will minimize the quasi-particle 
energy in the vortex state. 

Concluding Remarks. - We have analyzed the upper critical field Hc2 for the proposed 
p+h-wave superconductor PrOs4Sbi2 [6] and compared with the observed Hc2 (i) and H* it) . 
The present model reproduces Hc2 without any adjustable parameters. Also the Hc2(t) ob- 
served in the single phase crystals indicates the single phase is most likely the A phase. Also 
H*(t) appears to belong to the A phase. Does this suggest inhomogeneity? Also the present 
analysis predicts that the square vortex lattice is favored when the magnetic field is parallel to 
one of the crystal axes. So far no strong correlation between the superconducting transition 
temperature TcA and Tcb a,nd RRR [8] has been observed. This could be interpreted as 
meaning that disorder has little effect. However, from the early analysis of chiral p-wave su- 
perconductivity, we conclude that Hc2iT,T)/ Hc2iO,T) is rather universal. An open question 
is whether a high-quality single crystal with RRR > 100 would show a single or double phase 
transition. Clearly such a study is highly desirable. 
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